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Abstract 

Using 1.8 million DD pairs and a neutrino reconstruction technique, we have studied the de- 
cays D° -► K-e + v e , D° -► TT~e + u e , D + -► K°e + v e , and D+ -► 7r°e + ^ e . We find £>(D° -► 
Tr-e+^e) = 0.299(11)(9)%, £(£>+ -»• 7r°e+i/ e ) = 0.373(22)(13)%, B(D° -► K- e +i/ e ) = 3.56(3)(9)%, 
and B(D + — > K°e + z^ e ) = 8.53(13) (23)%. In addition, form factors are studied through fits to 
the partial branching fractions obtained in five q 2 ranges. By combining our results with recent 
unquenched lattice calculations, we obtain \V cd \ = 0.217(9)(4)(23) and \V CS \ = 1.015(10)(11)(106). 
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Study of the semileptonic decays of D mesons plays an important role in our understand- 
ing of the CKM matrix [lj]. These decays allow determination of \V CS \ and \V c d\ by com- 
bining measured branching fractions with results from unquenched lattice QCD (LQCD) 
calculations [2]. With improved branching fraction precision these measurements also pro- 
vide rigorous tests of LQCD Q. The tests can be approached by assuming unitarity of 
the CKM matrix and comparing the constrained matrix elements |4| to measured elements, 
or by comparing the measured and calculated ratios of semileptonic and purely leptonic 
branching fractions, which are independent of CKM matrix elements. Verification of LQCD 
calculations at the few percent level will provide validation for their application to the B 
system. 

This Letter presents a study of the D° — > K~e + u e , D° — > 7i^e + u e , D + — > K°e + i> e , and 
D + — ► 7r°e + z/ e decay modes (charge conjugate modes implied). A more detailed description 
of the analysis is provided in a companion article |5j. The results are based on 281 pb _1 
of e + e~ data at the ^(3770) resonance (1.8 million DD pairs) collected with the CLEO-c 
detector at the Cornell Electron Storage Ring (CESR) 0. The sample is a superset of the 
data used to obtain the first CLEO-c semileptonic branching fraction measurements |7| ■ For 
each mode we determine the partial branching fractions in five q 2 = m\ v ranges. Summing 
the rates yields the total branching, fraction; fitting the rates constrains form factor shapes; 
comparing to LQCD calculations [2j determines the CKM elements \V c d\ and \V CS \. 

The analysis technique rests upon association of the missing energy and momentum in an 
event with the neutrino four- momentum [8J], enabled by the hermeticity and excellent reso- 
lution of the CLEO-c detector. Charged particles are detected over 93% of the solid angle by 
two wire tracking chambers within a 1.0 T solenoid magnet. The momentum resolution is 
0.6% at 800 MeV/c. Specific ionization and a ring imaging Cerenkov detector (RICH) pro- 
vide particle identification; a CsI(Tl) electromagnetic calorimeter provides photon detection 
over 93% of Am and a tt° mass resolution of ~6 MeV/c 2 . 

Electron candidates are identified above 200 MeV/c over 90% of the solid angle by com- 
bining information from specific ionization with calorimetric, RICH and tracking measure- 
ments. To reduce our sensitivity to final state radiation (FSR), we add photons within 3.5° 
of the electron flight direction to the electron momentum. 

Signal charged pions and kaons are identified using specific ionization and RICH mea- 
surements. A 7T° candidate must have a 77 mass within 2.5 standard deviations (cr) of the 
7T° mass. Kg candidates are reconstructed using a vertex fit to candidate 7r + 7r~ daughter 
tracks. The 7r + 7r~ mass must be within 4.5a of the Kg mass. 

The missing four-momentum in an event is given by p miss = {E miss ,p miss ) = p total - 
SPcharged — DPneutrai, where the event four- momentum p to tai is known from the energy and 
crossing angle of the CESR beams. The charged and neutral particles included in the sums 
pass selection criteria designed to achieve the best possible |p m i SS | resolution by balancing 
the efficiency for detecting true particles against the rejection of false ones 

Association of the missing four-momentum with the neutrino is only valid if the event 
contains no more than one neutrino and all true particles are detected. We thus exclude 
events that have either more than one electron or non-zero net charge. The core |p m i SS | 
resolution in signal Monte Carlo (MC) simulated events satisfying these criteria is ~15 
MeV/c. 

To select signal events we require that the M^ iss = E^ iss — |p m i SS | 2 be consistent with 
a massless neutrino. Using MC simulations that included all resolution effects, candidate 
selection was found to be optimized by requiring |M^ iss /2|p miss | | < 0.2 GeV/c 3 . Since |p m i SS | 
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resolution is roughly half that of E miss , in subsequent calculations we take p v = ( |p m i SS | , Pmiss) • 
Semileptonic decays D — > /iez/ e , where h = tt or K, are identified using four-momentum 
conservation. Specifically, reconstructed D candidates are selected based on AE = (Eh + 
E e + Ey) — -Ebeam (expected to be close to zero within our 20 MeV resolution). The D 
momentum constraint is recast as the beam-constrained mass M\, c , defined below, which 
peaks near the D mass for signal. 

Selection criteria were optimized with independent MC samples. Background sources 
include events with hadrons misidentified as electrons (fake electrons), non-charm continuum 
production, and DD processes other than signal. The optimal energy requirement was 
determined to be —0.06 < AE < 0.10 GeV. For the Cabibbo-favored modes, the background 
remaining after this selection is only a few percent of the signal. For the Cabibbo-suppressed 
modes, significant background remains from signal-mode cross-feed and from the related 
modes D + — ► K^e + v e and D + — ► Kg(ir Tr )e + i> e , where the (7r°7r°) indicates the Kg decay 
mode. Restricting the AE of the non-signal side of the event reduces these backgrounds. We 
also require D + — ► 7r°e + z/ e candidates to have the smallest \AE\ of any final state candidate 
in the event, and that these events contain no reconstructed D° — > K~e + v e candidate. The 
average background fraction (g 2 -dependent) in the pion modes is about 20%. 

Since the \p v \ resolution dominates AE resolution, we improve our p v measurement by 
scaling it by the factor ( satisfying AE = (Eh + E e + (E v ) — E^, CSum = 0. We use (p u to 



calculate M bc = J E 2 cam — \p h + p e + (Pu\ 2 with a resolution of 4 MeV/c 2 . We calculate 
q 2 = (p v + p e ) 2 with a resolution of 0.01 GeV 2 /c 4 , independent of q 2 . 

To extract branching fraction information we perform a simultaneous maximum likelihood 
fit [o| to the Mb c distributions of the four signal modes in five q 2 ranges. The simultaneity 
automatically provides self-consistent rates for misreconstruction of one signal process as 
another (cross-feed) and for background from the related K° processes. The M^ c distribution 
is divided into 14 uniform bins over the range 1.794 < M^ c < 1.878 GeV/c 2 . To simplify 
the statistical interpretation of our results we limit the number of multiple entries per event: 
a given event can contribute to at most one D° and one D + final state. For multiple D° 
or D + candidates with M\ x > 1.794 GeV/c 2 , we choose the one with the smallest |Ai?|, 
independent of q 2 . 

We fit the data to the signal and five background components. The signal mode MC 
components are based on EvtGen 1(| with modified pole-model (BK parameterization) form 
factors 11] and parameters from the most recent unquenched LQCD calculation [2|. Several 
corrections, relating to inclusive D decay and reconstruction (see Ref. [!]), are applied to 
our GEANT-based [12j MC samples. These lead to few percent (or less) changes in the 
measured yields, and are determined precisely enough (using a large DD sample with one 
fully-reconstructed D meson per event) to yield sub-percent systematic uncertainties. We 
are also sensitive to the signal efficiency and kinematic distortions due to FSR. Based on the 
angular and energy distributions for FSR photons, we correct our signal MC, generated with 
PHOTOS [13|] without interference, to the leading-order KLOR [1^ calculations applied to 
charm decay. 

To reduce our sensitivity to form factors, we extract an independent rate for each of the 
five q 2 intervals corresponding to the reconstructed q 2 ranges in each mode (a total of 20 
yields) . 

We also use MC samples to describe the DD background and continuum contributions. 
We absolutely scale the continuum components according to their cross sections at the 
■^(3770) and the measured data luminosity. The non-signal DD sample was generated using 
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TABLE I: Branching fractions, isospin ratios, and form factor parameters (isospin corrected for the 
7T°e + f mode). Errors are (stat.)(syst.) or (stat.)(syst.)(theor.). Correlation coefficients (from com- 
bined statistical and systematic uncertainties) for variables in any two (three) preceding columns 
are given by p (pij). For the a« we have assumed |V^ S | = 0.976 and \V c d\ = 0.224. 



q 2 (GeV 2 /c 4 ) 

B(n-e+u e ){%) 

B(n°e+v e ){%) 

B(K'e+u e )(%) 

B(K°e+v e )(%) 

Ik 



< 0.4 0.4-0.8 0.8 - 1.2 1.2- 1.6 > 1.6 Total 

0.070(5)(3) 0.059(5)(2) 0.060(5)(2) 0.044(4)(2) 0.066(5)(2) 0.299(H)(9) 

0.084(10)(4) 0.097(H)(4) 0.062(9)(3) 0.063(10)(2) 0.067(H)(3) 0.373(22)(13) 

1.441(21)(35) 1.048(18)(28) 0.681(15)(18) 0.340(11)(10) 0.048(5)(12) 3.557(33)(90) 

3.436(82)(93) 2.544(73)(69) 1.589(58)(44) 0.821(42)(24) 0.139(18)(5) 8.53(13)(23) 

2.12(31)(9) 1.54(22)(7) 2.47(43)(13) 1.78(32)(7) 2.48(45)(13) 2.03(14)(8) 

1.06(3)(3) 1.04(4)(3) 1.09(5)(3) 1.05(6)(4) 0.88(15)(3) 1.06(2)(3) 



\V cq \ 

0.218(11)(5)(23) 
0.216(17) (6) (23) 
1.023(13)(13)(107) 
1.004(20) (15) (105) 









Series Parameterization 




Decay 


ao 


Ol 


02 


P01 


P02 P12 


|Vo,|/+(0) l + l/p-6 


P 


ir~e + u e 


0.044(2)(1) 


-0.18(7)(2) 


-0.03(35)(12) 


0.81 


0.71 0.96 


0.140(7)(3) 1.30(37)(12) 


-0.85 




0.044(3)(1) 


-0.23(H)(2) 


-0.60(57)(15) 


0.80 


0.67 0.95 


0.138(H)(4) 1.58(60)(13) 


-0.86 


K~e+u e 


0.0234(3)(3) 


-0.009(21)(7) 


0.52(28)(6) 


0.62 


0.56 0.96 


0.747(9)(9) 0.62(13)(4) 


-0.62 


K°e+v e 


0.0224(4)(3) 


0.009(32) (7) 


0.76(42)(8) 


0.72 


0.64 0.96 


0.733(14)(11) 0.51(20)(4) 


-0.72 



Decay 



K°e+u e 



|Vc,l/+(0) 
0.146(4)(2) 
0.149(6) (3) 
0.735(5)(9) 
0.710(8)(10) 



Simple Pole Model 
m P oic (GcV/c 2 ) p 
1.87(3)(1) 0.63 
1.97(7)(2) 0.65 
1.97(3)(1) 0.36 
1.96(4)(2) 0.53 



|Vc 9 |/+(0) 
0.142(4)(2) 
0.147(7) (4) 
0.732(B)(9) 
0.708(9)(10) 



Modified Pole Model 
a p 
0.37(8) (3) -0.75 
0.14(16) (4) -0.75 
0.21(5)(3) -0.42 
0.22(8)(3) -0.59 



EvtGen, with decay parameters updated to reflect our best knowledge of D meson decays. 
This component floats separately for each reconstructed final state, but is fixed over the 
five q 2 regions within that state, thereby reducing our sensitivity to inaccuracies in the 
D decay model. Finally, we input signal MC components for D + — > K^e + v e and D + — > 
Kg('K 7t )e + i> e , whose rates in each q 2 region are proportional to those for the reconstructed 
D+ -> K° s e + v e mode with K° s -> n+n-. 

The contributions of events with fake electrons are evaluated by weighting hadron- 
momentum spectra in candidate events with misidentification probabilities measured in other 
CLEO-c data. This component is included with a fixed normalization in the fit. 

We allow the fit to adjust the M DC resolution in the D° — > 7r~e + z/ e , D° — > K~e + u e , and 
D + 7r°e + z/ e modes by applying a Gaussian smear to these distributions. The result is 
that the signal MC M DC resolution in these modes, ~3.5 MeV/c 2 , is increased to match the 
data resolution of ~4 MeV/c 2 . 

The resulting M DC distributions, integrated over q 2 , are shown in Fig. [T]with the fit results 
(from D and D decays) overlaid. The value of the likelihood for this fit is — 2m£ = 275.5 
for 280 — 27 = 253 degrees of freedom. 

We obtain branching fractions (see Tabled]) for each q 2 region by combining the efficiency- 
corrected yields from the fit with the number of D°D° (N D om) and D + D~ {IS d+d-) pairs for 
our sample. An independent study of hadronic D decays [15J finds N d odo = 1.031(16) x 10 6 
and N D + D - = 0.819(13) x 10 6 . We find ratios of branching fractions R = B(D° — > -K~e + v e )/ 
B{D° -> K~e + v e ) = 8.41(32)(13)% and R+ = B{D + -> 7T°e + u e )/B(D + -> K°e + u e ) = 
4.37(27)(12)%. Tabled also lists the partial-width ratios 4 = T(D° -> n-e + u e )/T(D + -> 
7r°e + z/ e ) and I K = T{D° -> K-e + u e )/T{D + -> K°e + u e ), with lifetimes input from Ref. fl. 
Isospin symmetry predicts 1^ = 1 and Ik — 1- 

The systematic uncertainty (see Ref. 0) is dominated by uncertainty in the number of 
DD pairs and in neutrino reconstruction simulation. The latter includes inaccuracies in the 
detector simulation and in the decay model of the non-signal D. Mainly through use of 
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events with a reconstructed hadronic decay, we evaluate systematic bias, as a function of 
q 2 , for the efficiency of finding and identifying signal hadrons, identifying signal electrons, 
and fake electron rates. Similarly, uncertainties that affect the cross-feed rates, such as 
those associated with non-signal 7r° and tt~ production spectra, as well as K~ faking ir~, are 
also assessed as a function of q 2 . We correct statistically significant biases, and propagate 
the uncertainty of each study into our measurement uncertainty. The remaining systematic 
uncertainties include M^ c resolution, the effect of the single-electron requirement, MC FSR 
modeling, dependence on form factors, and the N DD determinations. 

Our primary form factor shape analysis utilizes a series expansion that has been widely 
advocated as a physical description of heavy meson form factors 17, 18, 19, 2p| : 
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The expansion results from an analytic continuation of the form factor into the complex 
t = q 2 plane, with a branch cut on the real axis for t > (Md + ^k^) 2 that is mapped by 
z(t, to) = (y/t + — t — y/t + — t )/(y/t + — t + y/t + — t ) onto the unit circle. The constants 
t± = (Md ± rriK,Tr) 2 , and t is the (arbitrary) q 2 value that maps to z — 0. The physical 
region is restricted to \z\ < 1, so good convergence is expected. P(q 2 ) accommodates sub- 
threshold resonances: P(q 2 ) = 1 for D — > 7r and P{q 2 ) = z(q 2 ) M D ,) for D — > K. The 
function <p(q 2 ,to) can be any analytic function. We report parameters that correspond to 



to = and the "standard" choice for (j> (see, e.g. Ref. |20j) that arises naturally in studies 
of unitarity bounds on J2 a t- 

For comparison purposes, we provide results based on the simple and modified pole mod- 
els [ll|. These parameterizations can typically accommodate the form factor shapes observed 
in previous measurements, but only with parameters that deviate from the underlying phys- 
ical motivation j2l|. Note that differing experimental sensitivities across phase space can 
result in differing parameter values for a non-physical parameterization. 

Each parameterization is fit to our measured rates for the five q 2 regions; parameter sys- 
tematic uncertainties are obtained from fits to the rates obtained for each systematic varia- 
tion. Tabled summarizes the results, and Fig. [2] compares fits for the three parameterizations 
in our most precise mode D° — > K~e + v e . For the series expansion, we also express our re- 
sults as physical observables: the intercept |V^ g |/ + (0) and 1 + 1//3 — <5oc df + /dq 2 \ q 2 =0 [20], 
which represents the effects of gluon hard-scattering (5) and scaling violations ((3). The 
results from D° and D + decay agree well. 

For the series expansion, our kaon data prefer a non-zero quadratic z term. The proba- 
bility of x 2 improves from 29% (22%) to 89% (44%) with that additional term for the K~ 
(K°) fit. The pion measurements lack the sensitivity to probe this term, and two and three 
parameter fits yield similar results for the first two parameters. Since a quadratic term 
appears preferred for the kaons, however, we include that term in our series fits to the pion 
data to improve the probability that our shape uncertainties bracket the true form factor 
shape. While three of the central values for 02 are an order of magnitude larger than the 
other terms, we stress that regions of parameter space with a 2 of similar magnitude to a 
and a\ fall well within the 90% hypercontour for the fit, so no conclusion can be drawn 
about the size of a 2 or (potential lack of) convergence of the series from these data. 

In the simple pole model, we fit for the intercept and the pole mass m po i e . In the modified 
pole model, we fix the leading pole mass at the physical value, and fit for the intercept and the 
parameter a, which determines the effective higher pole contribution. We obtain reasonable 
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FIG. 1: The reconstructed M^ c distributions, integrated over q 2 , for data (points) and components 
(histograms) from the fit (see text): signal MC (clear), cross- feed and non-signal DD MC (gray), 
continuum MC (light gray). The e + fake component (black) is negligible on this scale. 
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FIG. 2: Fits to the D° — ► K e + u e partial branching fractions using the three form factor models 
described in the text. All results are normalized to the series expansion result. 



X 2 values, but obtain pole masses that deviate from Mp* (Mjj*) in the kaon (pion) modes 
by over 3er for the most precise fits. The 1 + 1/(3 — 5 value for the series expansion fit to 
the K~e + v e data is over 3cr from the value of ~2 necessary for physical validity of the BK 
parameterization, while our values for the BK a parameters from the kaon modes imply 
1 + 1/(3 — 5 values tens of o away. 

We extract |V^| and |V^ S | by combining our |V^ g |/ + (0) results from the series expan- 
sion fit with the unquenched LQCD results g ff~**\o) = 0.64(3) (6) and f^ K) (0) = 
0.73(3) (7) (Table H]). Averaging the D° and D + results (heeding correlations), we find 
\V cd \ = 0.217(9)(4)(23) and (V^l = 1.015(10)(11)(106), with the /+ uncertainty listed last. 
The discretization uncertainty in the FNAL LQCD charm quark action dominates. 
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We also extract the ratio |V^d|/|V^ s | from the ratio of our measured form factors. Av- 
eraging over D° and D + modes, with correlations accounted for, gives \V C d\f^ D ~* n '{0)/ 
\V C s\f {D ^ K) (0) = 0.188(8)(2). Comparison to a recent light cone sum rules (LCSR) cal- 
culation where /^^(O) / f( D ~*^ (Q) = 0.84(4), yields \V cd \/\V cs \ = 0.223(10)(3)(11). 

In summary, we have measured branching fractions and their ratios for four semileptonic 
D decay modes in five q 2 bins. The branching fraction results are the most precise ever 
measured and agree well with world averages. Our modified pole a parameter results agree 



within 1.3a with previous determinations by CLEO III [23], FOCUS [24j, and Kev results 



from Belle [25|, but show over 3<r disagreement with Belle K\xv results and LQCD fits. The 
a parameters obtained with our individual Kev results are separated from the recent BaBar 
result [26[ by about 2.5cr. Our z expansion results agree with BaBar's at about the 2cr level, 
depending on the total level of correlation between the BaBar r\ and T2 parameters. We have 
made the most precise CKM determinations from D semileptonic decays to date, and the 
results agree well with neutrino-based determinations of \Vd\ and charm-tagged W decay 
measurements of \V CS \ jij. Overall, these measurements represent a marked improvement in 
our knowledge of D semileptonic decay. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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